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We examine the density profiles of the quantum degenerate Bose-Fermi mixture of 174 Yb- 173 Yb, 
experimental observed recently, in the mean field regime. In this mixture there is a possibility of 
tuning the Bose-Bose and Bose-Fermi interactions simultaneously using two well separated optical 
Feshbach resonances, and it is a good candidate to explore phase separation in Bose-Fermi mixtures. 
Depending on the Bose-Bose scattering length (Ibb , as the Bose-Fermi interaction is tuned the density 
of the fermions is pinched or swapping with bosons occurs. 
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I. INTRODUCTION 

The quantum degeneracy in a Boson-Fermi mixture 
was first experimentally realized for the system consist- 
ing of 7 Li and 6 Li [1] [2] . Since then it has been observed 
in several other Bose-Fermi mixtures, 23 Na- 6 Li [3], 87 Rb- 
40 K g], 87 Rb- 6 Li 0, 4 He- 3 He[5], 174 Yb- 173 Yb|3, and 
84 Sr- 87 Sr [5]. These are candidate systems to explore the 
effects of boson-induced fermionic interactions, of par- 
ticular interest is the boson mediated fermionic super- 
fluidity. Another property of interest is the dynamical 
instabilities of the fermionic component arising from the 
attractive fermionic interactions, which is also boson me- 
diated [TT]. Precondition to observe either of these is a 
precise control of the inter species interaction through a 
Feshbach resonance. Which has been observed in 87 Rb- 
40 K [9j [10] and used to trigger the dynamical collapse of 
40 K [TT] , the same is numerically analyzed in ref. [T2~] [J3] . 

Density distributions in the phase separated domain 
of the Bose-Fermi mixture is also an important property 
of interest. Like in binary condensates, dynamical insta- 
bilities can be initiated in the phase separated domain 
through manipulations of interaction strengths. For bi- 
nary condensates, the recent work on mixtures of two 
different hyperfine states of 87 Rb [14] is a fine example 
of controlled experiment on phase separation. 

In this regard, Molmer and collaborator [HI [16] ex- 
amined the zero temperature equilibrium density distri- 
butions and predicted widely varying density patterns 
as a function of inter species interactions. However, 
the bose-bose and bose-fermi interactions considered are 
extremely strong for experimental realizations. Simi- 
lar studies have examined the ground state geometry 
in spherical traps |171 118) . The conditions for mixing- 
demixing have been analyzed for homogeneous Bose- 
Fermi mixture [191 and Bose-Fermi mixtures inside traps 
[T71 |2"D] . Although, very high interaction strengths are 
achievable through magnetic Feshbach resonances in al- 
kali atoms, simultaneous tuning of both the boson-boson 
and boson-fermion interactions is not possible. How- 



ever, simultaneous tuning is possible with optical Fes- 
hbach resonances (OFR) when the resonant frequencies 
of the boson-boson and boson-fermion interactions are 
well separated. 

With the realization of quantum degeneracy in the 
mixture of 174 Yb- 173 Yb where intra-species interactions 
for 174 Yb can be tuned by OFR [2~T], we find it perti- 
nent to revisit these studies. With the possibility of tun- 
ing inter-species interactions for 174 Yb- 173 Yb mixture, 
it may be possible to realize the ground state geome- 
tries which are hitherto elusive. We, therefore, consider 
173 Yb- Yb mixture to study the density profiles for var- 
ious values of coupling strengths in the present work. It 
must also be mentioned that the isotopes of Yb exhibit 
a wide range of inter- and intra-species interactions. It 
has attracted lot of attention as selected isotopic com- 
positions may exhibit dynamical instabilities triggered 
through the interactions. The 174 Yb- 176 Yb is one such 
Bose-Bose binary mixture currently investigated for in- 
stabilities on account of the attractive intra-species in- 
teraction of 176 Yb [221 [23] . 

The paper is organized into four sections. In the next 
section, Section. II, we provide a brief description of the 
mean field equations of bosons and fermions. This is fol- 
lowed with the section on phase separation, where the 
nature of bose-fermi phase separation is discussed as a 
function of the inter-species interaction. More impor- 
tantly, the occurance of fermion pinching is explored. 
Position swapping between bosons and fermions, as the 
inter-species interaction is increased, is then examined in 
the next section. We then conclude with Section. V. 



II. ZERO TEMPERATURE MEAN FIELD 
DESCRIPTION 

We examine the stationary state properties of a Bose- 
Fermi (BF) mixture consisting of 174 Yb and 173 Yb in 
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spherically symmetric trapping potentials 



Vfe(r) = im B wV, Vp(r) 
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where the subscripts B and F stand for boson and 
fermion, respectively, and u is the radial trap frequency 
for the two components. The fermions are spin-polarized 
( single species ) and the fermion-fermion mean field 
interactions arise from the degeneracy pressure |24j . 
Whereas, the boson-boson and boson-fermion interac- 
tions arise from the s-wave scattering between the atoms. 
Considering these, the mean field energy functional of the 
Bose-Fcrmi mixture is [231 
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where ^(r, i) and ^ F (r,t) are bosonic and fermionic 
wave functions satisfying the normalization condition 
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Here Gbb = 47r/i 2 aBB^VB/TOB, where obb is the bosonic 
s-wave scattering length and Nb is the number of 
bosons, is the bosonic intra-species interaction, Gbf = 
2irh 2 aB F Np/mR and Gfb — ^f^dBF^B/mR, where 
mfl = (mBmp)/(mB + rn F ) is the reduced mass, N F 
is the number of fermions, and czbf is the inter-species 
scattering length, are the inter-species interactions, and 
A = h 2 (6n 2 N F ) 2 / 3 1 (2m F ) . The Lagrangian of the sys- 
tem 
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Using the action principle 
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we get a set of coupled partial differential equations 
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It is more convenient to rewrite Eqs. ^ in a di- 
mensionless form by defining dimensionless parameters 
in terms of the frequency u and the oscillator length 
aho = -\/^/( to bw). Using f = r/ctho, t — tw as the scaled 
dimensionless variables of length and time, respectively, 
Eqs. ^ can be rewritten as 
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where the rescaled wave functions are ipB — a t/^^B(i,i) 



and if) F — 
strength parameters are 
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Similarly, the interaction 
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with m rat io = m F /mB- For the sake of simplicity, we 
represent the scaled quantities without the tilde ( ~ ) in 
the rest of the article. For spherically symmetric trap- 
ping potential, the Eqs ^ are reduced to one dimensional 
coupled mean field equations 
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These are the coupled mean field equations which de- 
scribe the bose-fermi mixture in trapping potentials at 
zero temperature. To obtain the stationary solutions, we 
solve the equations numerically using Crank-Nicholson 
scheme |26j with imaginary time propagation. 



III. PHASE SEPARATION 

Broadly speaking, for large values of Gbb> the den- 
sity profiles of the boson-fermion mixture in spherical 
symmetric traps can have three distinct geometries in 
phase separated regime: (a) fermionic core surrounded by 
bosonic shell, (b) bosonic core surrounded by fermionic 
shell, and (c) shell of bosons between fermionic core 
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and fermionic outer shell [TS]. The inter-species inter- 
actions begin to play an important role, in determining 
the stationary state structure, when the density profiles 
of the bosons and fermions are of similar spatial extent. 
This occurs when the bosonic intra-species interaction 
is strong. In the Thomas-Fermi (TF) approximation, 
the necessary condition for a mixture of equal number 
of bosons and fermions (N B = N-p = N ) is [T5] 
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Here after, we use a BB to represent the particular value of 
obb at which bosons and fermions have the same spatial 
extent. For species which are isotopes of the same element 
the mass difference is small and 77i rat i ~ 1. The condi- 
tion is then reduced to 
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Considering N ~ I0 6 , which is the typical value in ex- 
perimental realizations, a BB ss 0.17ah o - As a^ is in 
general ~ 10 _6 m for harmonic trapping potentials, the 
required value of <z B b is in the strongly interacting do- 
main. It could be achieved when the interaction is tuned 
through a Feshbach resonance, magnetic in the case of 
alkali-Earth metal atoms. With the overlapping den- 
sity profiles, more intricate density patterns are observed 
when the inter-species interaction is increased, however, 
tuning q,bf with magnetic Feshbach resonance is ruled 
out. This complication does not arise when the inter- 
actions are tuned with well separated OFRs. In which 
case, the isotopes of two-valence lanthanide atom Yb is 
a suitable candidate. It has seven stable isotopes: five 
bosons ( 168 Yb, 170 Yb, 172 Yb, 174 Yb, and 176 Yb) and two 
fermions( 171 Yb, and 173 Yb ), and homo-nuclear OFRs of 
bosonic isotopes ( 172 Yb and 176 Yb ) were recently stud- 
ied [27]. Among the various possible species pairings 
174 Yb- 173 Yb, which has positive intra- and inter-species 
background scattering lengths, is an ideal candidate to 
study bose-fermi mixtures in the strongly interacting do- 
main. 

For our studies, we consider a 174 Yb- 173 Yb mix- 
ture containing 10 6 atoms of each species and trapped 
by spherically symmetric trap with trapping frequency 
w/(27r) = 400Hz. The obb is chosen to be equal 
to llOOao, which is achievable with the OFR of the 
6s 2 So — > 6s6p 3 P\ inter-combination transition. And, 
the bose-fermi inter-species scattering length obf can be 
tuned with OFR of the allowed 6s 2 1 Sq — > 6s6p 1 P\ tran- 
sition. This is a broad line and the disadvantage of using 
it is high atom loss rate. However, a major advantage 
of OFR tuned interactions is the fine spatial control it 
provides. Recently, submicron modulation of scattering 
length using Feshback resonances was achieved in 174 Yb 
|21j . Such precise control on the spatial variation of in- 
teraction strength is unrealistic with magnetic Feshbach 
resonances. From here on, where it is not explicitly men- 
tioned, reference to bose-fermi mixture implies 174 Yb- 
173 Yb isotope mixture. To examine the density profiles of 



the mixture in the strongly interacting domain, we keep 
obb fixed and vary obf so that the system progresses 
from mixing to full demixing regime via partial demixing 
regime. It must be emphasized that with TF approxima- 
tion, from Eq. ( 10 ) the spatial extent of density profiles 



with asF = are same when a BB = 1191. 71a . However, 
we have chosen a BB = H00a o , approximately the value 
at which the density profiles begin to exhibit the features 
of interest. 



A. Mixing to partial demixing regime 

Starting from the initial conditions of the mixture, 
which as mentioned earlier is equal spatial extent of the 
component species and obf ~ 0, the value of gibf is in- 
creased. To analyze the evolution of density profiles, con- 
sider the TF profile of the boson and fermions in scaled 
units as defined earlier 
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where uxy = 3xy /Ny and Ep is the fermi energy. From 
these expressions, the densities at the origin in absence 
of inter-species interaction are 
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1. Wave-function profiles of 174 Yb (dashed-green line) 
173 Yb (solid-orange line) in a mixture of 174 Yb- 173 Yb 
with Y B = N F = 10 6 , w/(2tt) = 400Hz and a BB = 1100a as 
cjbf is changed from mixing to partial demixing domain, (a) 
For a B F = 0.0, (b) for a B F = 0.0a BB and a BB = 1197. lloo, 
(c) for a B p = 0.6a BB , and (d) for a B F = 0.7a BB 



In TF approximation, the fermi energy and chemical 
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potential of the two species are 
E F =(6N) 1/3 , 
/i =i (l5a BB N) 2/5 . 
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Using these in Eq.(12|, the ratio of the densities at the 
origin is 
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Consider obb ~ 0.17, the value at which the profiles of 
the two species match for N = 10 6 . The population ratio 
at the origin is then 
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That is, when a BF = the fermion density is higher than 
the boson density at the center of the trap. Given this 
as the initial condition, when the inter-species interac- 
tion is switched on, the inter-species mean field energy 
G B f|*f(0)| 2 > Gfb|^b(0)| 2 . Hence, it is energetically 
favorable to shift the bosons from the center towards the 
edge of the trap. This is evident in the numerically ob- 
tained density profiles shown in Fig. [lja) , where the den- 
sity profile of the bosons is flattened around the origin 
and has higher density at the edges. 
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FIG. 2. Interspecies interaction energy E^t between 174 Yb 
(boson) and 173 Yb (fermion) as a function of (Jbf. The max- 
ima of -Bint occurs at <xbf = 0.7cibb. 

The inter-species interaction energy in mixing regime 



enters the partial demixing regime and the characteristic 
signature of which is a maxima in inter-species interac- 
tion energy. For our present calculations, the variation 
of the inter-species interaction energy with obf is shown 
in Fig. [2j The condition for attaining partial demixing 
in spherical traps is [20] 
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The above condition for partial demixing is evaluated us- 
ing TF approximation for the density profiles of both the 
components. From equation Eq.(16), the critical value 



of a BF required to reach partial demixing regime for the 
Bose-Fermi mixture under consideration is 0.44ob and 
is significantly lower than the value of 0.7a BB obtained 
from numerical solution of the coupled mean field equa- 
tions Eq([8]). The difference may be attributed to the 
simplifying assumptions in deriving the location of the 
inter-species interaction energy extrema. One of which is 
choosing the density profiles at a BF — to calculate the 
inter-species interaction energy. At a BF — 0.7obb there 
is dramatic decrease in the density of the Bosons near 
the trap center. This is accompanied by corresponding 
decrease in overlap region between the two components 
as is shown in Figjljb). 
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3. The wave function profiles of 174 Yb (dashed-green 
and 



173 Yb (solid-orange line) in a mixture of 174 Yb- 
N F = 10 6 , a B s = HOOao and oj/(2tt) = 



line) 

173 Yb, with iV B 
400Hz. For (a) obf = 0.75obb and for (b) obf = l.OaBB- 
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where we have used n B m n F « A^/(47ri? 3 r F/3) [2Q] with 
i?TF = \/2(6A r ) 1 / 3 for the system considered in the 
present work. As obf is increased further the system 



B. Partial demixing to phase separation 

A further increase of obf, beyond the critical value, 
enhances the segregation of the two species. This lowers 
the inter-species overlap and balances the larger interac- 
tion energy from higher obf ■ Ultimately, at higher values 
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of obf the overlap is almost zero, the system can then 
be considered fully phase separated. The condition to 
attain phase separation or fully demixed regime is [2"0] 
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For 174 Yb- Yb mixture with the previously mentioned 
parameters, the above criterion translates into obf > 
0.9aBB- In the phase separated domain, the separation 
occurs around the inner point where densities are equal. 
To identify the location of this point consider the obf = 
density profiles. If the two profiles intersect at r\, then 
from the TF approximation r t is the solution of the equa- 
tion 
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For the system of our interest Vp and Vb are al- 
most identical. Further, when obb is chosen (it satisfies 
Eq. (10)) to match the spatial extents of the densities, 
Ep ~ /i. Following which, to a very good approximation 
[E F - V F (n)} « \n- V B (n)]. The solution of Eq. ^TJ is 
then 
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The importance of r; is for the following: np(r) > 
n B {r) for r < n, and np-(r) < ns(r) for r > r\. For 
the 174 Yb- 173 Yb mixture, based on the above relation 
n = 13.09a ho for N B = Np = 10 6 and a BB = 0.17a ho , 
while the numerical value is r\ — 12.42ai 10 - Energetically, 
when obf is switched on it is favorable to accommodate 
the bosons and fermions at the outer and inner regions 
about n, respectively. As <xbf is increased, the position 
of r; tend to migrate outward but not dramatically. 

With further increase in obf, the bosons are expelled 
towards the edge of trapping potential while fermions 
are squeezed towards the trap center (see Figg. With 
TF approximation, the effective potential experienced by 
173 Yb in the overlap region is 
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1 for 174 Yb- 173 Yb 
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where we have considered m ra tio 
mixture. Obviously, the effective potential experienced 
by fermions vanishes at obf — clbb 1 an d this explains the 
constant wave function profile of 173 Yb in the overlap 
region as is shown in Fig[3](b). Unlike two component 
BECs, the density of the fermions in the region occupied 
by bosons is not zero when the criterion for full demixing 
is satisfied. 
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4. The wave function profiles of 174 Yb (dashed-green 
line) and 173 Yb (solid-orange line) in a mixture of 174 Yb- 
173 Yb, with iVe = N-p = 10 6 , a BB = 1100a and cj/(2tt) = 
400Hz, as obf is steadily increased from an initial value of 
<jbf = Abb - For (a) obf = 1.05aBB, for (b) obf = 1.1(Jbb, for 
(c) obf = 1.15aBB, and for (d) obf = 1.25aBB- 



C. Fermion pinching 



For the values of Obb marginally below Og B , besides r; 
there is another point r Q where the densities are identical. 
The location of r a is rather sensitive to kinetic energy 
corrections of the bosons 129 



■Ekir 



2.5§LlB 



(24) 



Without the kinetic energy correction, that is with TF 
approximation, r exists up to higher values of obf- How- 
ever, the kinetic energy correction softens the profile at 
the edges and r vanishes as q,bb approach a BB . In the 
phase separated domain when r is close to the edge, the 
fermion density is depleted at r\ < r < r Q for higher obf- 
And, there is fermion density enhancement at r < r\ and 
r > r a . For the bosons it is opposite: there is density 
enhancement at rj < r < r a , and depletion at r < r\ and 
r > r a . 

As <zbf is increased to values larger than o,bb, the ef- 
fective potential within the overlap region (r; < r < r G ) 
is approximately 
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where 77 = |m ra ti — 3bf/5bb| and like in the previous case 
we can take m rat i ~ 1. The form of V g is repulsive with 
a maxima at r; and decreases towards r Q . The net effect 
is, the fermion density profile is pinched at the region 
where r is marginally larger than r\. Onset of pinching 
is clearly discernible in Fig. |4^a) and Fig. [4jb) shows 
density profile with even higher pinching effect. At higher 
values of obf the pinching is complete and an island of 
fermions appears at the edge. The Fig. |4jc-d) show the 
formation of the fermionic island due to pinching in the 
173 Yb- 174 Yb mixture considered in the present work. 
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IV. PROFILE SWAPPING 

A remarkable feature in the evolution of density pro- 
files as a function of obf is the observation of profile 
swapping for certain range of parameters. In which 
the fermions are initially at the core and bosons form 
a shell. However, at higher values of obf the bosons 
occupy the core and fermions forms a shell around it. 
As an example to illustrate profile swapping, consider 



xi(r 




15 o 5 

r (a OS c) 

FIG. 5. The wave function profiles of 174 Yb (dashed-green 
line) and 173 Yb (solid-orange line) in a mixture of 174 Yb- 
173 Yb, with Nb — Nf — 10 5 , a B s = 1100a and u/(2n) = 
400Hz, as obf is steadily increased from an initial value of 
<2bf = ass- For (a) (Jbf = LOobb, for (b) obf = 1.05aBB, for 
(c) Obf = l.laBB, and for (d) a B F = 1.15obb- 



N B = N F = 10 5 , from Eq. (10) the spatial extents 
are equal at gibb = 0.24aho- However, retain the value 
<2bb = H00a o as in the case of 10 6 atoms in each 
species. In this case, the spatial extent of the bosons 
is less than the fermions, however, there are two points 
at which the densities of the bosons and fermions are the 
same. As mentioned earlier the ground state geometries 
of Bose-Fermi mixtures in spherical symmetric traps can 
be broadly categorized in three types, and this is evi- 
dent from the Figj6]for 174 Yb- 173 Yb mixture with with 
-/Vb = N-p = 10 6 and w/(2tt) = 400Hz. 

When <2bf is set to a non-zero value, at lower values 
the changes in the equilibrium density profiles exhibit 
a pattern similar to fermion pinching. Like in fermion 
pinching as <zbf is ramped up, there is a depletion of 
fermions from the overlap region as shown in Fig. [5ja- 
d). However, at some value of obf a dramatic departure 
occurs. The fermions from the core are expelled to the 
edges and bosons settle at the core, Fig [7] At interme- 
diate values of obf, the bosons form a shell sandwich 
between fermions at the core and an outer shell. This is 
evident from the density profiles shown in Fig. [5jci) . As 
is evident from the figures, the migration of the fermions 
to the flanks occurs at a relatively minute change in obf, 
from 1.15aBB to 1.16obb- 

To analyze the profile swapping based on total energy 
considerations, take the density profiles just prior to the 
expulsion of fermions from the core. The inter-species 
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FIG. 6. The phase diagram of 174 Yb- 173 Yb mixture with 
Nb = Nf — 10 6 and uj/(2tt) = 400Hz. Dashed-green and 
solid-orange curves are the semi-analytic conditions for mix- 
ing to partial demixing and demixing to phase separation (or 
full demixing) transitions, while dotted-blue is the numeri- 
cally obtained criterion for mixing to partial demixing transi- 
tion. In phase separated regime, for asB < 750 bosonic core 
is surrounded by fermionic shell, for 750 < obb < 1100 there 
is fermion pinching along with profile swapping at obb ~ 750. 
For 750 < cibb <■ 1100, bosonic shell is surrounded by 



fermionic core and shell. For qbb 
is surrounded by bosonic shell. 
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interaction energy is 
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Here, r x like defined earlier is the inner point where the 
two densities are equal. The first term is the inter-species 
interaction energy arising from the inner boundary of 
the overlap region. And S is the interpenetration depth 
considered symmetric for simplicity. The second term is 
the interaction energy from the remaining overlap region, 
though the upper limit of integration is taken as oo in re- 
ality it extends up to the point where tib is nonzero. To 
simplify the analysis assume that the fermions from the 
core, after the position swapping, are pushed beyond the 
overlap domain. The interaction energy when swapping 
occurs is 



(irUBF^B^-F 



n+S 



druBF^B^F, (27) 



where in the first term, the lower limit accounts for the 
nonzero fermion density around the core. The occurancc 
of position swapping implies that 



dr"UBF?lB«F > / drUBF^B"-F, (28) 

n-S JO 



at some value of obf- In other words, like in binary 
mixtures of condensates, at some point the geometry of 
the overlap region determines the nature of the density 
profile. 
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FIG. 7. The wave function profiles of 174 Yb (dashed-green 
line) and 173 Yb (solid-orange line) in a mixture of 174 Yb- 
173 Yb, with A^b = Nf ~ 10 5 , a BB = 1100a and u/(2tt) = 
400Hz, as (Zbf is steadily increased from its initial value of 
<2bf = 1.16aBB- For (a) obf = 1.16aBB, for (b) «bp = 
1.17aBB, for (c) cebf = 1.18aBB, and for (d) obf = 1.2cibb. 



Profile swapping, initiated by tuning interspecies scat- 
tering length, appears to be a promising tool to study 
Rayleigh- Taylor type of instability in Bose-Fermi mix- 
tures, which has been theoretically studied in two species 
Bose-Einstein [SUl GH] condensates. The idea is to start 
with a ground state geometry with fermions forming the 
core, and then increase obf so that the new ground state 
has the fermionic core swapped by bosonic one. 



V. CONCLUSIONS 

We have analyzed the equilibrium density profiles of 
the 174 Yb- 173 Yb Bose-Fermi mixture for a range of in- 
teraction strengths. In this Bose-Fermi mixture, it is 
possible to tune both the Bose-Bose and Bose-Fermi in- 
teractions across a range of values. Density profiles of 
the two species display pinching and position swapping 
when the boson-boson scattering length is close to a BB , 
the value at which the spatial extent of the bosons is 
same as the fermions. Pinching occurs when the obb is 
marginally below a BB , at these values as g,bf is increased 
fermions at the edges are pinched to form a thin shell, 
whereas at even lower values of obb , as asp is increased 
the fermions are expelled to the edge and density profiles 
are swapped. At intermediate values of a B p the profiles 
undergo through a series of configurations, and these are 
significantly different from the ones in Bose-Bose mix- 
tures. Close to the profile swapping domain, it should be 
possible to initiate Rayleigh- Taylor instability through a 
controlled variation of obf- This would be significantly 
different from Rayleigh- Taylor instability in condensates. 
In future, it would be interesting and important to ex- 
plore various instabilities which may occur at the Bose- 
Fermi interface boundaries. These could be qualitatively 
different from their analogues of binary condensates. 
ACKNOWLEDGMENTS 

We thank S. A. Silotri, B. K. Mani, and S. Chattopad- 
hyay for very useful discussions. The numerical computa- 
tions reported in the paper were done on the 3 TFLOPs 
cluster at PRL. 



[1] A. G. Truscott, K. E. Strecker, W. I. Mc Alexander, 
G. B. Partridge, R. G. Huletdagger, Science 291 2570 
(2001). 

[2] F. Schreck, L. Khaykovich, K. L. Corwin, G. Ferrari, 

T. Bourdel, J. Cubizolles, and C. Salomon, Phys. Rev. 

Lett. 87, 080403 (2001). 
[3] Z. Hadzibabic, C. A. Stan, K. Dieckmann, S. Gupta, 

M. W. Zwierlein, A. Gorlitz, and W. Ketterle, Phys. Rev. 

Lett. 88, 160401. 
[4] G. Roati, F. Riboli, G. Modugno, and M. Inguscio, Phys. 

Rev. Lett. 89, 150403 (2002). 
[5] C. Silber, S. Giinther, C. Marzok, B. Deh, 

Ph. W. Courteille, and C. Zimmermann, Phys. Rev. 

Lett. 95, 170408 (2005). 
[6] J. M. McNamara, T. Jeltes, A. S. Tychkov, W. Hoger- 

vorst, and W. Vassen, Phys. Rev. Lett. 97, 080404 

(2006). 

[7] T. Fukuhara, S. Sugawa, Y. Takasu, and Y. Takahashi, 

Phys. Rev. A 79, 021601(R) (2009). 
[8] M. K. Tey, S. Stellmer, R. Grimm and F. Schreck, Phys. 

Rev. A 82, 011608 (2010). 
[9] F. Fcrlaino, et ai, Phys. Rev. A 73, 040702(R), (2006); 

S. Inouye, J. Goldwin, M. L. Olsen, C. Ticknor, 

J. L. Bohn, and D. S. Jin, Phys. Rev. Lett. 93, 183201, 



(2004). 

[10] M. Zaccanti, C. DErrico, F. Ferlaino, G. Roati, M. In- 
guscio, and G. Modugno, Phys. Rev. A 74, 041605(R) 
(2006). 

[11] G. Modugno, G. Roati, F. Riboli, F. Ferlaino, 
R. J. Brecha, and M. Inguscio, Science 297, 2240 (2002). 

[12] M. Modugno, F. Ferlaino, F. Riboli, G. Roati, G. Mod- 
ugno, and M. Inguscio, Phys. Rev. A 68, 043626 (2003). 

[13] S. K. Adhikari, Phys. Rev. A 70, 043617 (2004). 

[14] S. Tojo, Y. Taguchi, Y. Masuyama, T. Hayashi, H. Saito, 
and T. Hirano, Phys. Rev. A 82, 033609 (2010). 

[15] K. Molmer, Phys. Rev. Lett. 80, 1804 (1998). 

[16] N. Nygaard and K. Molmer, Phys. Rev. A 59, 2974 

(1999) . 

[17] A. Minguzzi and M. P. Tosi, Phys. Lett. A 268, 142 

(2000) . 

[18] Robert Roth, Phys. Rev. A 66, 013614 (2002). 

[19] L. Viverit, C. J. Pethick, and H. Smith, Phys. Rev. A, 

61, 053605 (2000). 
[20] Z. Akdeniz, A. Minguzzi, P. Vignolo, and M. P. Tosi, 

Phys. Rev. A 66, 013620 (2002). 
[21] R. Yamazaki, S. Taie, S. Sugawa, and Y. Takahashi, 

Phys. Rev. Lett. 105, 050405 (2010). 



8 



[22] K. Kasamatsu and M. Tsubota, J. Low Temp. Phys. 150, 
599 (2008). 

[23] G. K. Chaudhary and R. Ramakumar, Phys. Rev. A 81, 
063603 (2010). 

[24] D. A. Butts and D. S. Rokhsar, Phys. Rev. A 55, 4346 
(1997). 

[25] P. Capuzzi, A. Minguzzi, and M. P. Tosi, Phys. Rev. A 

67, 053605 (2003). 
[26] P. Muruganandam and S. K. Adhikari, Comp. Phys. 

Comm. 180, 1888 (2009). 



[27] K. Enomoto, K. Kasa, M. Kitagawa, and Y. Takahashi, 

Phys. Rev. Lett. 101, 203201 (2008). 
[28] M. Kitagawa, et al, Phys. Rev. A 77, 012719 (2008). 
[29] F. Dalfovo, L. Pitaevskii, and S. Stringari, Phys. Rev. A 

54, 4213 (1996). 
[30] S. Gautam, D. Angom, Phys. Rev. A, 81, 053616 (2010). 
[31] K. Sasaki, N. Suzuki, D. Akamatsu, and H. Saito, Phys. 

Rev. A. 80, 063611 (2009). 



